lattice Boltzmann simulation of electrowetting-on-dielectric in a rough-wall channel 
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A lattice Boltzmann model was proposed to simulate electrowetting-on-dielectric (EWOD). The insulative 
vapor and the electrolyte liquid droplet were simulated by the lattice Boltzmann method respectively, and the 
linear property between cosine of contact angle and the electric field force confirms the reliability of this model. 
In the simulation of electrolyte flowing in a rough- wall channel under an external electric field, we found that a 
narrow channel is more sensitive than a broad channel and the flux decreases monotonously as the electric field 
increase, but may suddenly increase if the electric field is strong enough. 
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With the development and progress of nano-science and 
nano-technology, how to manipulate tiny amounts of liquids 
on surfaces becomes great challenge to nano-engineering. 
Gradients of both temperature and concentration of sur- 
factants across droplets may produce forces to propel 
droplets yj], but compared to these methods, chemical and 
topographical structuring of surfaces are much finer for con- 
trolling the equilibrium morphology of a droplet. To further 
actively control the wettability of a droplet, we may coat the 
droplet with self-assembled monolayers IID, but considering 
the switch speed, range of contact angle variation, and long 
term stability, electrowetting (EW) may be the best method to 
control a droplet on a surface. 



In electrowetting, electrowetting on dielectric (EWOD) [2] 
and electrowetting on insulator coated electrodes (EICE) yfl 
must be the two widely used methods to control a droplet, 
because of the several advantages: without any electrochem- 
ical reactions, clean surfaces ensure reversibility, and highly 
nonwetting surfaces may be used to enhance electrowetting 
effect [4]. Recently, the experimental research on electrowet- 
ting has made great progress: electrowetting has been pro- 
moted to aplication of lab-on-a-chip |j5J] ; through an external 
voltage, adjusting hydrophobic nature of an oil droplet con- 
taining dissolved dye, an electrowetting-based reflective dis- 
play has been invented [6]; changing the contact angle of a 
sessile droplet via electrowetting, shape of the droplet, as well 
as the focal length of this droplet can be changed, and a flex- 
ible liquid lenses may be made |7]. Theoretically, the con- 
tact angle should decrease rapidly as the external voltage in- 
creases, but in experiment, the contact angle has been found to 
be saturated when the external voltage is high enough [ 2]. Be- 
cause many factors may affect the contact angle of a droplet in 
experiment, the saturation mechanism of this angle keeps un- 
known lyj], despite several various mechanisms have been pro- 
posed [ 1] . In order to study the variation of contact angle on 



electrowetting as well as the electric effect on microchannel 
and microfluidics, it is desirable to propose a numeric method 
to study electrowetting of tiny electrolyte droplet on an insu- 
lative layer. 

Since the lattice Boltzmann method (LBM) has been 
proposed, it has succeeded in simulating complex fluid flows 
including multi-phase flow [fToh particle suspension flow (HI 
d, binary mixture HQ and blood flow QEl! H]. Be- 
cause of its intrinsically parallel and simpleness, it has been 
recognized as an alternate method for computational fluid dy- 
namics 11711 . Through introducing intermolecular forces, the 
lattice Boltzmann method may be used to s tudy microfluidics 
lfl8ll . microchannel and nanochannel flows [19]. 



In this Letter, base on the multiple phases lattice Boltz- 
mann model lllOh . considering the electric force exerted on the 
electrolyte near an insulative layer, a lattice Boltzmann model 
was proposed to investigate two-dimensional electrolyte elec- 
trowetting on a smooth insulative layer and the effect of 
EWOD on the electrolyte flow in a rough- wall channel. 

In two dimensions, the system of the D2Q9 lattice Boltz- 
mann model Jsh is described by a single particle distribu- 
tion function / a (x, t), which can be thought as the number 
of fluid particles at site x, time t, and possessing one of the 
nine velocities e a with eo = (0,0), e a = (cos7r(a — l)/2, 
shi7r(a — l)/2),a = 1,2,3,4, and e a = \/2(cos7r(2a — 
l)/4,shi7r(2a - l)/4), for a = 5,6,7,8. The distribution 
function evolves according to a Boltzmann equation that is 
discrete in both space and time, 



/ a (x + e a , t + 1) - / a (x, t) = —(f a - f e J). (1) 

r 

The density p and macroscopic velocity u are defined as 

P = ^2fa, pu = ^2f a e a . (2) 
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In order to reproduce the Navier-Stokes equation, the equi- 
librium distribution function f^ q may have the following 
form li 
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f? = w a p[l + 3e a • u + ^ (e a • u) 2 - (3) 

where wq = 4/9, = 1/9, for a = 1, 2, 3, 4, and w a = 
1/36, for a = 5, 6, 7, 8. 

Considering the intermolecular interaction, through an 
interaction potential among particles in the LBM, Shan 
and Chen improved the LBM to simulate multiple phases 
fluid OH: 

V(X,X') = G a (X,X , )^(x)^(x / ), (4) 

where G a (X, X') is a Green function and G a (X, X') = Q for 
|X - X'| = y/2, G a (X,X') = 2Q for |X - X'| = 1 and 
G CK (X,X / ) = otherwise. Here we take <p(x) = po(l — 
exp(—p/ po)) with po of a constant. The resulting force is 

F CT (x) = -(p(x) ^2 G a (p(x + e a )e a . (5) 

If the node at (x + e a ) is a solid node, similar to a fluid node, 
the according resulting force can be written as 112011 

¥ w (x) = -if(x) W a S(x + e a )e a , (6) 

a 

where W a = w, for |e a | = a/2, W a = 2w for |e a | = 1 and 
W a = for |e a | = 0. If the node at (x + e a ) is a solid node, 
5 = 1, else 5 = 0. 

For an EWOD system, which consists of a conductive 
droplet and a metallic electrode insulated with an insulative 
layer, the double conductive layer has a fixed capacitance per 
unit area, c = SoSd jd.d and Sd are the thickness and dielectric 
constant of the insulative layer, eo is the vacuum permittivity. 
The free energy of the system reads 

F e = A iv ai v + A sv cr sv + A s i(a si - °* ) - ApV", (7) 

where A is the area of the interface and a is the interfacial ten- 
sion. The subscript Iv, sv, and si indicate the liquid-vapour, 
solid- vapour, and solid-liquid interfaces respectively. U is the 
potential difference between the conductive droplet and the 
metallic electrode. Ap describe the pressure drop across the 
liquid- vapour interface and V is the volume of the droplet. 
Through minimizingthe free energy F e , we obtain the basic 
equation of EWOD W 

cos(^)-cos(^o) = ^^, (8) 
zaaiv 

where 8 is the contact angle of the droplet and #o is that 
without external electric feild. The electric force exerted on 
the bottom layer of the conductive droplet per unit area is 
f = fn = £Q£ 2d n an d n is the unit normal vector of the 
layer pointing to the metallic electrode. Thus equation ^ is 
simplified as 

cos(0) - cos((9 ) = — . (9) 



Based on equation ©, through a force exerted on the single 
layer of the liquid near the insulative layer, we may simulate 
an electrolyte droplet electrowetting on an insulative layer by 
the LBM. 

Fig. [T]is a droplet adopted on a smooth substrate simulated 
by the LBM. The single layer of fluid near the substrate is 
shown in Fig. [2] The fluids in the center with a high density 
pi consist the liquid, the other fluids in the two sides with a 
low density p v belong to the vapor. There is a transition part 
connect the liquid and the vapor. Because the liquid is con- 
ductive but the vapor is insulative, it is necessary to introduce 
a critical density p c (p v < p c < pi) to distinguish the fluid 
from conductor (with p > p c ) and insulator (with p < p c ). In 
Fig. [2 the fluids in the range between x\ and X2 in the sin- 
gle layer near the substrate are conductive. Fluids outside this 
range are regarded as insulative vapor. 

Under the influence of such forces, the momentum at each 
fluid node is changed to 

pu / = pu + r(F a +F w +f), (10) 

where u' is the new velocity used in equation (Q}. 




FIG. 1 : Shape of a droplet without external electric field. 

We now present the simulation results on a box consisting 
of 100 x 50 lattice units, we fixed the fluid-fluid interaction 
strength Q = -0.52, fluid- wall interaction strength w = —0.14 
and the relaxation time r = 0.9. The constant po in the ex- 
pression of potential ip was chosen to be 1 .0. At beginning, the 
fluid was homogenous and possessed a density of 0.85. An at- 
tractive force was assumed to push the fluids to condense into 
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FIG. 2: The density of the single layer of fluid near the substrate, pi 
and p v are the density of the conductive liquid and insulative vapor 
respectively. p c is a critical density to distinguish whether the fluid is 
conductive (p > p c ) or insulative (p < p c ). The fluid in the range of 
xi < x < X2 is regarded as conductive. Here we choice p c = 1.0. 
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a droplet in the middle of the substrate. When the droplet 
had formed, the artificial attractive force was eliminated and 
a normal droplet formed in the middle of the substrate shown 
as in Fig. [TJ Next, an electrostatic force of attraction / was 
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FIG. 3: cos(#) — cos(#o) versus electric force /. and 0o are the 
contact angle of the droplet with and without external electric field. 



ulate EWOD. Fig. [3] shows the contact angle changes with the 
electric force. From this figure we can see that cos 6 — cos #o 
increases proportionally with electric force /, consistent with 
Eq. ©. Thus, this simple model can be used to simulate 
EWOD. 





FIG. 4: Sketch of rough-wall channel. 




11111111 





FIG. 5: Rough- wall channel with width w — 160. The body force 
pointing right is f b = 8.0 x 10~ 6 . The electric force: (a) / = 0.0, 
(b)/ = 0.15, (c)/ = 0.25. 

exerted on the single layer of fluids near the substrate to sim- 



FIG. 6: Rough- wall channel with width w = 80. The body force 
pointing right is f b = 8.0 x 10" 6 . The electric force: (a) / = 0.0, 
(b)/ = 0.15, (c) / = 0.25. 
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FIG. 7: Average of fluid velocity along x-axis u x in the rough- wall 
channel with width w = 80. 



For a nanochannel, as the channel size decreases, the 
surface-to- volume ratio increases. Therefore, various prop- 



erties of the walls, such as surface charges and roughness, 
greatly affect the fluid motions in the nanochannels [21] . Re- 
cently, suffice roughness effect on micro-flow has been widely 
studied from theory | 22 |] . molecular dynamic (MP) and other 
numeric methods 112 ll l23l 12411 and experiment 112511 . Here, 
we focus on the effect of EW on an electrolyte flow inside 
a rough- wall channel. Fig. [4] shows a sketch of rough- wall 
channel. The top and the bottom walls of the channel were 
both equipped with steps of width a and height b. The cen- 
treline distance between two neighbouring steps was p. The 
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FIG. 8: Relative flux of the channel versus electric force. Q and 
Qo are the electrolyte flux in the channel with and without external 
electric field. The rightward body force exerted on the electrolyte is 
f b = 8.0 x 1(T 6 . 

channel width was w and the left and right boundary of the 
channel were periodic. The electrolyte was insulated from 
the channel. Applying a voltage between the channel and the 
electrolyte, we can used this system to investigate EW effect 
on the electrolyte flow in a rough- wall channel. Fig. [5] shows 
the electrolyte flowing through the rough- wall channel under 
a rightward body force /& with different external electric field. 
From this diagram we can see that when the applied voltage 
increases, the electric force / increases (/ ex U 2 ), the elec- 
trolyte wets the wall, and the electrocappirary phenomenon 
is obviously seen. For a more narrow rough-wall channel, as 
shown in Fig. [6l when the electric force is strong enough, the 
balance of the system will be break down and the electrolyte 
can not stay in the center of the channel. Fig. 0(c) shows that 
the electrolyte are absorbed into one rough-wall of the chan- 
nel living only little parts of the electrolyte keep contacting 
with another rough- wall. Although the electrolyte keeps flow- 
ing rightward, the contact points keep static, because of small 
velocity of the flow. 

For a narrow channel, not only the properties of the wall 
affect the flow very much, to an electrolyte flow, but also the 
external electric field or charges. Fig. [7] shows the average of 
fluid velocity along x-axis u x under different external electric 
fields. The according dimensionless flux Q/Qo is shown in 
Fig. [U Q is defined as: 

Q = i J u x dxdy, (11) 

where L is length of the channel and Qo is the flux without 
extenal electric field. From this figure we can see that the 
dimensionless flux in a narrow rough- wall channel decreases 
more rapidly than in a broad channel. When the electric field 
is strong enough, for a narrow rough-wall channel, the elec- 
trolyte leans towards one wall and lives thick vapor layer be- 
tween the electrolyte and another wall (see Fig. [6] (b)) and 
the vapor layer makes the wall slip more (shown in Fig. [7] 



for / = 0.25). Thus the flux increases suddenly under such 
external electric field (see Fig. [S]for / = 0.25). 

In conclusion, we introduced a critical density p c to distin- 
guish conductive liquid from insulative vapor and an electric 
force exerted on the conductive liquid to improve the LBM 
to simulate EWOD. The linear property of contact angle con- 
firmed the reliability of this moel. Further, we applied this 
model to simulate electrolyte flowing in a rough-wall chan- 
nel and found that the external electric field affect the density 
distribution, velocity and flux greatly for a narrow rough- wall 
channel. For a narrow rough-wall channel, when the external 
electric field is strong enough, the electrolyte may be adopted 
into one rough- wall and departure another rough- wall, reduc- 
ing the whole resistance and resulting a suddenly increasing 
of flux in the channel. 
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